Binary black hole coalescence in the extreme-mass-ratio limit: 
testing and improving the effective-one-body multipolar waveform 
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We discuss the properties of the eflective-one-body (EOB) multipolar gravitational waveform 
emitted by nonspinning black-hole binaries of masses and M in the extreme-mass-ratio limit 
^i/M — V <^ 1. We focus on the transition from quasicircular inspiral to plunge, merger and 
ringdown. We compare the EOB waveform to a Regge- Wheeler- Zerilli waveform computed using 
the hyperboloidal layer method and extracted at null infinity. Because the EOB waveform keeps 
track analytically of most phase differences in the early inspiral, we do not allow for any arbitrary 
time or phase shift between the waveforms. The dynamics of the particle, common to both wave- 
generation formalisms, is driven by a leading-order 0{v) analytically resummed radiation reaction. 
The EOB and the Regge- Wheeler-Zerilli waveforms have an initial dephasing of about 5 x 10~* rad 
and maintain then a remarkably accurate phase coherence during the long inspiral (~ 33 orbits), 
accumulating only about —2 x 10"'^ rad until the last stable orbit, i.e. A0/0 ~ —5.95 x 10~®. We 
obtain such accuracy without calibrating the analytically-resummed EOB waveform to numerical 
data, which indicates the aptitude of the EOB waveform for studies concerning the Laser Interfer- 
ometer Space Antenna. We then improve the behavior of the EOB waveform around merger by 
introducing and tuning next-to-quasicircular corrections in both the gravitational wave amplitude 
and phase. For each multipole we tune only four next-to-quasicircular parameters by requiring com- 
patibility between EOB and Regge- Wheeler-Zerilli waveforms at the light ring. The resulting phase 
difference around the merger time is as small as ±0.015 rad, with a fractional amplitude agreement 
of 2.5%. This suggest that next-to-quasicircular corrections to the phase can be a useful ingredient 
in comparisons between EOB and numerical-relativity waveforms. 

PACS numbers; 04.30.Db, 95.30.Sf, 



I. INTRODUCTION 

In the last few years numerical and analytical relativ- 
ity have demonstrated how to use information from the 
strong-field-fast-motion regime of coalescing black-hole 
binaries to build accurate analytical models of their dy- 
namics and of the gravitational radiation emitted . 
Although numerical- relativity (NR) simulations of binary 
black holes have reached a high degree of accuracy and 
flexibility p^ - [l6l |. a comprehensive spanning of the mul- 
tidimensional parameter space remains prohibitive. Ana- 
lytical models are thus of fundamental importance to set 
up the bank of gravitational wave (GW) templates for de- 
tection. The limiting case is given by extreme-mass-ratio 
inspirals (EMRIs) and mergers; NR simulations simply 
can not access such regime and post-Newtonian (PN) 
techniques are inaccurate at such velocities. We need an- 
alytical models for the GW emission from EMRI systems 
because they are primary target sources for the Laser In- 
terferometer Space Antenna (LISA) and because their 
parameter space is very large [17|. 

The only analytical approach currently capable of ac- 
curately following the complete dynamics and providing 
waveforms (inspiral-plunge-merger-ringdown) of coalesc- 
ing black-hole binaries is the effective-one-body (EOB) 
approach to the general relativistic two-body dynam- 
ics [ll|i 118142^ . The EOB formalism employs resummed 
PN results (for dynamics and waveforms) in order to ex- 



tend their validity in the strong-field-fast-motion regime, 
i.e. in a region where they are inaccurate in their standard 
Taylor-expanded form. In brief the analytical construc- 
tion is based on (i) a dynamics governed by a resummed 
Hamiltonian and an expression for the mechanical an- 
gular momentum loss (the radiation reaction) and (ii) 
a waveform-generating algorithm which combines a pre- 
scription to resum the Taylor-expanded PN multipolar 
waveform up to the merger and a matching procedure to 
the quasinormal-mode (QNM) waveform to describe the 
postmerger phase (an oscillating black hole). 

One key aspect of the EOB approach is its flexibil- 
ity [2^. Although the formalism is based on analytical 
results known only at a given PN order, it is possible 
to take into account (yet uncalculated) higher-order ef- 
fects by means of suitable flexibility parameters. These 
parameters may be determined (or just constrained) by 
comparison with results from numerical-relativity simula- 
tions valid in the strong-field-fast-motion regime. Several 
recent works 0, H, 0-11, [l^l have shown how this tuning 
can be implemented to obtain analytical waveforms that 
match the numerical ones within numerical errors. The 
tuned EOB formalism can then be used for parametric 
studies. 

The Regge- Wheeler- Zerilli (RWZ) metric perturbation 
theory is the natural tool to compute the GW 

emission from a system of two nonspinning black holes, 
of masses /i and M, in the extrcme-mass-ratio limit 
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(EMRL) V = /i/M ^ 1. In this regime several numerical 
results can be used to calibrate the EOB dynamics and 
waveforms [sl. [28l - [3]| . In particular, recent gravitational- 
self-force calculations [12, helped in putting con- 
straints on the functions entering the EOB conservative 
dynamics 0, [s^l- The Regge- Wheeler- Zerilli perturba- 
tion theory has been used for many years in the Fourier 
domain (see, for example, (ssj and references therein) and 
neglecting radiation-reaction effects, since Davis, RufHni, 
and Tiomno computed the waveform emitted by a par- 
ticle radially plunging into the black hole [s^. Only 
recently the RWZ approach has been extensively devel- 
oped in the time rfomam I37H42I with the inclusion of the 
radiation-reaction force |32l. l33l. lisl - lisj . 

Time-domain simulations using the perturbation the- 
ory are efficient and accurate and complement NR sim- 
ulations in EMRL. The first calculation of the com- 
plete gravitational waveform emitted during the transi- 
tion from inspiral to plunge, merger and ringdown in the 
EMRL was performed in Ref. [43[, thanks to the combi- 
nation of the RWZ perturbation theory and the 2.5PN 
accurate (analytical) Pade-resummed radiation-reaction 
force [23. Reference [2^ used that result as a target 
waveform to assess the performances of the correspond- 
ing EOB (resummed) analytical waveform. The compar- 
ison was restricted to the quadrupole case, m = £ ^ 2. 
The knowledge from that study was useful in subse- 
quent EOB/NR waveform comparisons. The treatment 
of the analytical radiation reaction in the strong-field- 
fast-motion regime has been improved since then, thanks 
to a resummed and factorized form of the PN multipolar 
waveform [lllilii]. 

In [ill (hereafter paper I) two of us presented an ac- 
curate computation of the gravitational radiation gen- 
erated by the coalescence of two circularized nonspin- 
ning black holes in the EMRL. The results were ob- 
tained with an improved version of the finite-difference 
code of [2^ [4^ , which implements the expression of the 
radiation-reaction force based on the (5PN-accurate) an- 
alytical waveform resummation of [^ . The knowledge of 
the "exact" RWZ multipolar waveform opened the way 
to two main conclusions, extensively discussed in paper I: 
first, the computation of the final kick velocity imparted 



to the system by GW emission, v /{cv ) 



0.0446. 



This value proved consistent with the corresponding one 
extrapolated from a sample of numerical-relativity sim- 
ulations [l^l (see Fig. 7 and Tables IV and V in paper 
I), as well as with the outcome of an independent cal- 
culation that relies on a different treatment of the ra- 
diation reaction [4J]. Second, it was possible to show 
a very good agreement (at the 10~^ level) between the 
mechanical angular momentum loss provided by the an- 
alytical expression of the radiation reaction and the GW 
angular momentum flux computed from the RWZ wave- 
forms. This second result supports the consistency of 
our approach. Notably, the agreement between the two 
functions was excellent also below the last stable orbit 
(LSO) and almost along the entire plunge phase up to 



merger (see Figs. 8 and 9 of paper I). The results of [45| 
also turned out to be compatible with the first NR com- 
putation of binary black hole coalescence in the large- 
mass-ratio regime (1:100) [l^. Recently [i^ we further 
improved the RWZ approach of paper I by combining it 
with the hyperboloidal layer method [i^ . This approach 
brings two main benefits. First, it allows us to extract 
GWs at null infinity thereby eliminating the gauge 

effects related to the GW extraction at a finite radius. 
In addition, because we evolve the RWZ equations on a 
smaller coordinate domain, we substantially improve the 
efficiency of our code. 

The aim of this paper is to perform a detailed compar- 
ison, multipole by multipole, between the RWZ and the 
corresponding analytical waveforms computed within the 
EOB approach. For the particle dynamics both codes 
(RWZ and EOB) iinplement the resummed radiation- 
reaction force J^^ of [31 updated to include 5PN-accurate 
terms also for subdominant multipoles. The latter come 
from the 5.5PN-accurate (Taylor-expanded) circularized 
multipolar waveform computed by Fujita and Iyer [46j . 
The particle dynamics is computed within this 5PN- 
accurate (resummed) approximation and is the same in 
both codes. For simplicity, we decided not to improve it 
further by tu ning the resummed flux entering the radia- 
tion reaction jSOl f3l| . For the waveform we compare the 
full multipolar structure up to m = ^ = 4, going beyond 
the simple quadrupole contribution. 

The waveform comparison brings new knowledge with 
respect to the flux comparison of paper I (see also [2§|) 
for two main reasons. First, we assess the performance 
of the resummed EOB waveform in describing the phase 
of each multipole. Second, we perform detailed analyses 
of the ncxt-to-quasicircular (NQC) corrections that are 
needed in the late-plunge phase. NQC effects are actu- 
ally responsible for the differences in the EOB and RWZ 
fluxes in the strong-field-fast-motion regime, as it was 
pointed out in paper I (see Fig. 8 there and also the re- 
lated discussion in 2^). At the waveform level, several 
studies [1, [3, [13, [Is I have demonstrated that NQC cor- 
rections to the EOB waveform (and radiation reaction) 
are needed to improve its agreement with the numerical 
one during the late-plunge and merger phase. Previous 
works were restricted to the quadrupole case. Two cen- 
tral new benefits of this paper are (i) the assessment of 
the complete multipolar EOB waveform in the EMRL 
during the transition from inspiral to plunge and merger; 
and (ii) the development of a robust procedure to tune 
NQC corrections to the gravitational wave amplitude and 
phase. 

The paper is organized as follows. In Sec.|ll]we summa- 
rize the main features of our RWZ numerical target wave- 
form described in detail elsewhere [45|, [i^. In Sec. lIIII we 
describe the structure of the multipolar EOB waveform 
of [^ [4^, giving all the details of the implementation 
used here. In Sec. IIVI we first present our results for the 
inspiral phase and then describe the procedure to tune 
NQC parameters necessary to improve the EOB wave- 
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form at merger. The discussion is based mainly on the 
£ = 2 multipoles. In Sec. |V] we assess the quality of the 
complete EOB multipolar waveform, discussing explic- 
itly multipoles up to ^ = 4. We finally put together some 
concluding remarks in Sec. I VII Two appendixes are in- 
cluded to complement the information given in the main 
text. Throughout this paper we use geometrized units 
with c = G = 1. 



II. REGGE-WHEELER-ZERILLI WAVEFORMS 




We compute numerical waveforms at future null in- 
finity via the time-domain RWZ pcrturbative approach 
introduced in [i^ and improved in [H, HI]. We per- 
form a hyperboloidal evolution of the RWZ equations 
with a point-particle source modeling the smaller-mass 
black hole. The distributional S function representing the 
particle is approximated by a narrow Gaussian of finite 
width a <^ M . The dynamics of the particle is started us- 
ing post-circular initial data as defined in [l^ , which 
generate negligible eccentricity at the beginning of the 
evolution. The conservative part of the dynamics is gov- 
erned by the — ?> limit of the EOB Hamiltonian (the 
Hamiltonian of a particle on Schwarzschild spacetime) 
with the following, dimensionless variables: the relative 
separation r = R/M, the orbital phase (p, the orbital 
angular momentum p^p — Pip/{fiM), and the orbital lin- 
ear momentum pr, = Pr, / ^J■, canonically conjugate to the 
tortoise radial coordinate separation r* = r+2 ln(r/2— 1). 
The expression for the analytical radiation-reaction force 
is described in @, and has been updated with 
the new terms in the resummed waveform amplitude re- 
cently computed in [i^ at fractional 5PN accuracy. In 
the computation of we sum over multipoles up to 
i = 8 included. The dynamics is then computed by solv- 
ing Eqs. (l)-(7) of paper I. 

Because of the different analytical approximation to 
the fiux and because of the hyperboloidal evolution, the 
RWZ waveforms employed here are quantitatively new 
with respect to those of paper I. By contrast, at a qual- 
itative level, there are no appreciable differences. In the 
following we shall comment only on the main new fea- 
tures of our pcrturbative approach. For a complete de- 
scription of the method, the equations, and the notation, 
we refer to jH El, IH . 

We adopt the hyperboloidal layer method for the RWZ 
equations to extract waves at future null infinity and to 
increase the efficiency of our code j49l - [5l| . The essential 
ingredient of this approach is a suitable transformation of 
the standard Schwarzschild time coordinate t in combi- 
nation with spatial compactification. We anticipate here 
part of the technical steps discussed in [i^ . 

We solve the RWZ equations, written in a general co- 
ordinate system [s^l , using coordinates (t, p) S R+ x 
The coordinates (r, p) coincide with {t,r^,) 
in a domain D~^'^ = [R~ , , that entirely includes the 
motion of the particle. The RWZ equations on -D~_^ have 




{t-S)/M 

FIG. 1: Multipolar "convergence" of the TZh+ / {AI i/) polariza- 
tion of the Regge-Wheeler-Zerilli multipolar waveform. Top 
panel: The complete wave train (~ 37 cycles). Bottom panel: 
Impact of subdominant multipoles around the merger time. 
The vertical dashed line indicates the light-ring crossing. 



then the same form as in paper I. The spatial compact- 
ification sets in at the interface p = = i?+ in a suffi- 
ciently differentiable way. The compactifying coordinate 
p(r*) maps the infinite domain = [i?+,oo) to the 
finite p domain = , S] , where S > Rf is a con- 
stant. A new time coordinate is introduced according 
to the prescription that the timelike Killing field is left 
invariant, dt = 9r, which implies 

T = t-/i(n). (1) 

The height function /i(r*) is related to p(r*) by the con- 
dition that the representation of outgoing null rays is left 
invariant 

t - = r - p . (2) 

Equations ^ and Q imply /i(r*) = r* — p{r^) that, 
together with the choice of a sufficiently smooth spatial 
compactification, determines a (future) hyperboloidal fo- 
liation. Thus the surface p = S corresponds to future 
null infinity I"*", and outgoing waves are evenly resolved 
in the compactifying coordinate p. Our numerical do- 
main reads [R~,S]j^+ = [—50,70] with the interface at 
i?+ = 50 and is covered by TV = 3001 grid points 

The RWZ-based approach we use to compute the nu- 
merical target waveforms relies on certain approxima- 
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tions. Our first approximation is to consider the dissi- 
pativc radiation reaction at leading order in the mass 
ratio ly, neglecting higher-order corrections that enter 
both the conservative and the nonconservative parts of 
the dynamics [H, [H, [11] . If one is interested in com- 
puting very long inspiral waveforms, these higher-order 
effects must be properly taken into account for LISA- 
rclated data analysis (for example, using the EOB for- 
malism [13, [m [Hi)- The computation of long inspiral 
waveforms with the 0{v) radiation reaction is affected by 
systematic uncertainties that depend on the mass ratio ly 
and on the integration time. By contrast, when one fo- 
cusses only on the late-time part of the waveform, i.e. the 
one corresponding to the "quasigeodesic" plunge, merger, 
and ringdown, one can extrapolate finite-z/ results to the 
V = limit. In paper I wc followed this procedure for 
dynamical quantities like the kick velocity or the energy 
emitted during the plunge, and wc showed that the im- 
pact of dealing with a finite value of v becomes less and 
less important when v < 10^^. Our second approxima- 
tion is that the expression for the radiation-reaction force 
J^^ is based on resummed PN results for circularized bi- 
naries and, as such, it neglects nonquasicircular correc- 
tions that explicitly depend on the radial momentum and 
its derivatives. The accuracy of this second approxima- 
tion can be checked only a posteriori but it is typically 
quite good also during the plunge [H, [i^ . 

On the basis of the considerations of paper I, we choose 
the mass ratio ly = 10~^ and the initial separation 
rg = 7M for the present EOB and RWZ comparison. 
These values guarantee a long inspiral (~ 37 orbits) as 
well as an accurate representation of the late-time wave- 
form (see its mild dependence on u when moving from 
v = 10~^ and 10~^ in Fig. 4 of paper I). We follow 
here the same normalization convention for the RWZ 
^£m°^ master functions [with (e)ven-parity and (o)dd- 
parity modes] as in paper I, i.e. 

(3) 

where TZ is the distance from the source, ^max is the 
maximum number of multipoles that we consider, and 
_2yfm = -2y£m(07*i') arc the s = —2 spin-weighted 
spherical harmonics computed in the convention of 
Ref. [13] ■ Figure [1] displays various multipolar approx- 
imations to the TZh+/{Mi/) polarization of the RWZ 
waveform along the fiducial direction (6,$) = (7r/4,0). 
The waveforms are shown versus retarded time at 
T — S. The dashed line in the figure (red online) refers to 
the complete waveform obtained summing the multipoles 
up to ^ = 8. The other two lines depict the partial con- 
tributions to the total waveform up to ^ = 4 (solid line) 
or £ = 6 (dash-dotted line) and neglecting the to = 
multipoles. The bottom panel of the figure is a close-up 
on the waveform around the conventional "merger" time, 
i.e. the time ^lr at which the particle crosses the light- 
ring r = 3M (vertical dash-dotted line) . The figure gives 



us a visual idea of the impact of the various multipoles on 
the accuracy of the total waveform, and it indicates that 
^max = 4 (with only m ^ 0) gives us a good approxima- 
tion, especially up to iLR- For this reason, and to keep 
the discussion sufficiently simple, we shall fix £niax = 4 
and consider only m 7^ multipoles to compute the RWZ 
target waveform to be compared with the EOB analytical 
waveform. Wc will discuss the fine details of the complete 
RWZ waveform in [i^ (see also [isj). 

III. EFFECTIVE-ONE-BODY RESUMMED 
MULTIPOLAR WAVEFORM 

In this section we review the structure of the EOB 
waveform [3, [1, [ll|, [4^. The EOB-resummcd mul- 
tipolar waveform can be split into two parts: the 
inspiral-plus-plunge-and-merger (insplumerg) waveform 
computed during the dynamics of the particle up to 
merger and the ringdown waveform that describes the 
post merger waveform as a superposition of black-hole 
QNMs. A simplified and efficient representation of the 
transition between the late-plunge and the ringdown 
regimes is accomplished by matching the insplumerg 
waveform to the ringdown waveform^. The complete 
EOB multipolar waveform reads 

hf^^ (^) {t; af™, a± ) = 0(t,„ - t) /^'-Pi™ (^) (t; a^) 

(4) 

where e denotes the parity, i.e. even (e = 0) for mass gen- 
erated multipoles and odd (e = 1) for current generated 
ones^. Since the particle motion is planar, e is equal to 
the parity of the sum (£-|-to). Furthermore, we explicitly 
highlight the dependence on the NQC parameters of™ 
and the QNM (complex) frequencies cr^^. 
The ringdown waveform is written as 

^ringdown (e) ^ ^ C^^^e'^^^ + ^ C" ^e"^-*, (5) 
n n 

where, following Ref. [1^, we use the notation a^^ = 
One ± iwn^ for the positive and negative QNM frequen- 
cies and C^^j^ for the corresponding amplitudes (note 
that, for simplicity, wc omitted here the parity index 



^ We introduced the nomenclature "insplumerg waveform" to in- 
dicate the part of the EOB waveform that is usually called in- 
splunge waveform in the literature [^. IllH . The reason for this 
choice is twofold: first, because it includes NQC corrections that 
become relevant essentially only around the merger time, and 
second, because it needs compatibility conditions with the RWZ 
waveform around merger to be fully determined. 

^ For notational consistency with previous analytical work 0| we 
label even- and odd-parity modes with e = and e = 1, re- 
spectively, while in Eq. l[3} above we used the labelling (e) and 
(o). 
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e). Here and ant indicate the frequency and the 
inverse damping time of each mode respectively, and 
ri = 0, 1, 2, . . . — 1 label the overtone number (n = 
denoting the fundamental mode). 

The insplumerg waveform can be written as the prod- 
uct of several factors. We factorize the NQC correction 
as 

^in^pWrg (e) ^ ^i„splu„go(.) ^NQC ) ^ (g) 

and the insplunge waveform is given as the product of the 
Newtonian contribution and a PN (resummed) correction 

by 



, insplunge (e) _ J, (^^0 



The Newtonian contribution is given by 



(7) 



where Y^™'{9,(t)) are the usual scalar spherical harmon- 
ics (computed on the equatorial plane 9 = tt/2), and 
the numerical coefHcients and ci+e(y) are explicitly 
given by Eqs. (5)- (7) of Ref. [l]. Because we work here 
in the EMRL [v — > 0), we pose Cf+e(0), leaving only 
the overall factor v. Following Ref. (consistently with 
Rcf . |2S . 1 5 Si] ) , the argument x in the Newtonian prefactor 
of Eq. dXl) is taken as 

x = vl = (rny , (9) 

where is the orbital frequency. This choice is preferable 

to X - — c i r c 

= ^circ '^'^^ to the violation of the circular 
Kepler's law during the plunge phase. The quantity 



(10) 



represents a factorizcd (and resummed) version of all the 
PN corrections to the waveform. It is given as the prod- 
uct of four factors: the /i-normalized effective source S^g , 
the tail factor r*"™ that resums an infinite number of 
leading logarithms entering the tail effects, the supple- 
mentary phase Sim, and the residual modulus correction 
P£m ■ The even-parity effective source S^g is given by the 
/i-normalized Hamiltonian of the system computed along 
the dynamics, while the odd-parity one S^g is given by 
the corresponding (Newton-normalized) angular momen- 
tum. The explicit expression of the tail factor as a func- 
tion of the orbital frequency reads 



T{£+1- 2ik) 
r(£+l) ' 



(11) 



where k = mfi, k = Mk. Following [4y|, we denote by ros 
the quantity that was previously denoted as tq in Q , and 
we choose tqs = = 2M / ^/e for consistency between the 
computations of the phase of [1^ and This value is 
chosen to match PN results (in harmonic coordinates) 



with black-hole perturbation results in Schwarzschild co- 
ordinates used here. 

Finally, the factor /ifo^^(af™) is a correction that mod- 
els the noncircular effects in the waveform. This effec- 
tive term allows us to analytically compute the wave- 
form beyond the circular approximation. The necessity 
of using a NQC corrective factor to the waveform am- 
plitude (usually denoted as f^'~^'^) was first pointed out 
in [29[ and then used (with variations) in several stud- 
ies [1 i, 11, S The inclusion of such a NQC cor- 
rection proved necessary to improve the closeness of the 
771 = £ = 2 FOB waveform to the numerical one during 
the late-plunge and merger dynamics. Building on this 
knowledge we extend the use of h^^'^ [1, to the other 
multipoles, with the important difference that we con- 
sider also corrections to the waveform phase 0,1111. This 
generalized NQC factor takes the form 



, NQC. fmN _ A , „lm Pr, 



X exp 



1 a; 



im Pr, 



(12) 



where the af"*'s are NQC flexibility parameters that have 
to be determined, multipole by multipole, by impos- 
ing some compatibility conditions between the FOB and 
RWZ waveforms around the merger time. 



IV. COMPARING EOB AND RWZ 
WAVEFORMS 

In this section we present the comparison between 
RWZ and FOB waveforms. In the first part we assess 
the performances of the resummed insplunge waveform, 
Fq. (O, to describe the inspiral up to the LSO cross- 
ing. The insplunge waveform includes neither NQC cor- 
rections nor matching to QNMs. Nonetheless it shows a 
very good agreement with RWZ during the whole inspiral 
up to (and even below) the LSO. In the second part we 
focus on the transition from inspiral to plunge, merger 
and ringdown, and we compare to the RWZ waveform 
the three different analytical representations of the re- 
summed multipolar waveform introduced above. Specif- 
ically we consider the same insplunge waveform men- 
tioned above; the insplumerg waveform that includes 
only the NQC corrections (to both phase and ampli- 
tude), Fq. ©; and the full FOB waveform, Fq. (g]), with 
NQC corrections and with QNMs matching. We em- 
phasize the necessity of introducing NQC corrections to 
the resummed FOB waveform and we propose compat- 
ibility conditions with the RWZ waveform to determine 
the NQC parameters of". The procedure we discuss is 
equally robust for all multipoles and does not require any 
relative time and phase shifting of the waveforms, or any 
hand-tuning of parameters. 

We present our comparison consistently with the no- 
tation of paper I; i.e. we use "Zerilli-normalized" metric 
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TABLE I: Phase diflFerence A(?!)'"=p'""«5c 
relative amplitude difference (AA/A)'"=p'"" 



EOB 

'Im 



j,RWZ 
J EOB 



and 



Afrn )/Aern between insplunge EOB and RWZ waveforms 
computed at to — 500M and at the LSO crossing, tLSO = 
4076. IM, for all the (m / 0) multipoles with 2 < ^ < 4. 



e m 




V^LSO 


?o 


\ / /500J\/ 




2 1 


-3.91 xlO 


-4 


-4.04x10" 


-3 


-1.16x10"'* 


-9.62x10"" 


2 2 


4.16 xlO' 


-4 


-2.48x10" 


-3 


-1.40x10"^ 


-2.72x10"=* 


3 1 


-1.74 xlO 


-3 


-1.86x10" 


~2 


-5.56x10"^ 


-4.6x10"^ 


3 2 


-3.29 xlO 


-4 


-4.65x10" 


-3 


7.54x10"* 


9.54x10"* 


3 3 


3.19 XlO" 


-4 


-3.76x10" 


-3 


-1.98x10"^ 


-4.10x10"^ 


4 1 


-2.06 xlO 


-3 


-2.36x10" 


-2 


-2.96x10"'* 


-8.02x10"^ 


4 2 


-1.35 XlO 


-3 


-1.51x10" 


-2 


-9.55x10"'' 


-3.98x10"^ 


4 3 


2.78 XlO" 


-3 


5.69x10^ 


4 


1.15x10"^ 


1.80x10"^ 


4 4 


4.22 XlO" 


-4 


-4.40x10" 


-3 


3.42x10"^ 


7.1x10"=* 



multipoles "^im- 
multipoles as 



These are related to the TZhim metric 



(13) 



where Ng = 1/ + 2){£ +!)£{£ - 1) and the label X 
stays for either insplunge, insplumerg, or EOB. Note that 
this convention implies a phase shift of tt/2 between the 



odd-parity multipoles 7^/l^2 and "ifYJ^. The EOB and 
RWZ multipolar waveforms are computed from the same 
dynamics; the RWZ one is typically extracted at and 
shown versus the corresponding retarded time t — S. The 
EOB waveforms is parametrized by the dynamical time 
t that will always be used as the reference time axis. 
Finally, the phases and the amplitudes Aim of the 
(RWZ or EOB) complex numbers are defined with 
the convention \I'p„, = A/ime""^^™ . 



(1) 



A. Quasiadiabatic inspiral 
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FIG. 2: Testing the waveform resummation for I = 2 at the 
beginning of the inspiral. Top panel; Insplunge and RWZ 
waveforms extracted at . Bottom panel: The phase dif- 
ferences A./.fO'^f^^^ = .^fo*^ - for RWZ waveforms 
measured at , are contrasted with the corresponding ones 
for RWZ waveforms measured at a finite extraction radius 
rf = lOOOM. 



Let us focus first on the quality of the resummed wave- 
form during the long (~ 37 orbits) quasiadiabatic inspi- 
ral. We identify here the end of the inspiral (and the 
beginning of the plunge) as the time tLSO = 4076. 27\/ at 
which the particle crosses the last stable orbit tlso = 
6M. This is clearly a convention because the transition 
from the inspiral to plun ge i s a blurred process that oc- 
curs around tlso uM, UM, l56[ . 

The upper panel of Fig. [5] displays the early-time evo- 
lution of the real part of the £ ^ 2 EOB multipoles 
(dashed lines) (top m — 1; bottom m — 2) together with 
the RWZ ones (solid lines). After the initial unphysi- 
cal transient the plot shows a remarkably good agree- 
ment between phases and amplitudes of the two wave- 
forms. The time evolution of the corresponding phase dif- 
ferences (during the complete inspiral) A<j>fl^^^'^'^ (t) = 



•^fm^l^) ~ '^^^(*) shown (as thicker lines, red onhne) 
in the lower panel of Fig. [2] The vertical dashed line 
in the lower panel of the figure (with label tLSo) indi- 
cates the LSO-crossing time. At the beginning, say for 
t/M < 500 (first 4 orbital cycles), the magnitude of both 
A(/)f2°^^^^(t) and A0fOBRWZ(^) ^elow lO'^. More 
precisely, at t = io = 500M, we have A(j)f^^^"^'^ {to) = 
4.16 X 10"4 and A(j)f^^^'^^ (to) = -3.98 x 10"! Such a 
small dcphasing is compatible with the expected uncer- 
tainty on the 4.5PN-accurate phases Sgm in Eq. pO|) . By 
evaluating the last (the 4.5PN ones) terms in Eqs. (5.8a)- 
(5.8b) of [401 at the initial position tq = 7M, we estimate 
54.5PN _ ^ X S4.5PN _ 5 27 X 10"^. Even if 

the phase difference grows by an order of magnitude on 
the interval [ioj^LSo]j it remains small at tLSO- Table |T] 
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accounts for the complete information, multipolc by mul- 
tipole, about phase and fractional amplitude differences, 
{/^A/A)ira = (AfOB _ ^RWZ)/^RWZ^ ^oth at to and at 

*LSO- 

When we sum together all the multipoles in the com- 
plete waveform, Eq. ([3]), we initially obtain A0gQQ'^'y"^° = 
5 X 10~^ (corresponding to total GW phase (j) = 51.74 
rad), that becomes as large as —2 x 10"'^ at the LSO 
{4> = 472.1 rad). This means that the complete EOB in- 
splunge waveform has dephased from the numerical one 
by only about -2.5 x 10"^ rad over the 420 rad (~ 33 
orbits) of GW phase evolution on the interval [toj^LSo], 
which yields A0EOBRWZ/0RWZ ^ .5 95 ^ iQ-e 

We point to one of the main findings of Ref. [i^ in the 
lower panel of Fig. [21 The figure depicts the phase dif- 
ferences to numerical waveforms extracted at the finite 



radius 



= lOOOM and displayed versus the corre- 



sponding observer retarded time. The plot clearly shows 
that the values of A(j)f^^^'^ obtained with waves ex- 
tracted at a finite radius (though very large) are approx- 
imately 2 orders of magnitude larger than those obtained 
with waves extracted at null infinity. The phenomenon 
discussed here in the £ = 2 case remains the same for 
the other multipoles. In particular, the largest and the 
smallest dephasing is for the m = I and the m = £ multi- 
poles respectively. The dephasing might grow (up to the 
~ 0.1 rad level) for > 2. The small dephasing of the 
waveform from the hyperboloidal layer calculation is the 
result of accurate wave extraction at null infinity (48| . 
At the present stage, this simple plot teaches us a very 
useful lesson: Extracting waves at lOOOM, even in the 
EMRL, may introduce dephasings that are considerably 
larger than the ones expected on the basis of the analyt- 
ical knowledge of the GW phase. A priori one expects 
such a phenomenon to be equally relevant, or even more 
important, in the comparable mass-ratio case due to the 
smaller PN accuracy at which the Sim^s are known when 
^ 0. 



B. Compatibility conditions at merger 

Let us focus now on the late-time part of the wave- 
form corresponding to plunge, merger and ringdown, 
i.e. t > ^Lso- The discussion mainly focusses again on 
the £ — 2 modes as a paradigmatic example. The rel- 
evant information is collected in Fig. |31 The top pan- 
els show the various waveform moduli divided by v, 
i.e. A2m/i^ = |^'2m|/i^7 with self-explanatory labelling. 
The bottom panels show the corresponding instantaneous 



^ We checked this result also using the numerical setup of paper 
I, i.e. with evolution along Cauchy time surfaces on a finite- 
size r* domain r« G [—500, 1500] and artificial boundaries (with 
Sommerfeld's outgoing boundary conditions) instead of a hyper- 
boloidal layer. A thorough discussion of those effects will be 
given in Ref. l48l . 



TABLE II: Strong-field-fast-motion information from the 
RWZ waveform at the light-ring crossing t — t^n = 4308. 4M 
and used to determine the af™ coefficients via the conditions 
given by Eqs. (fT4)l - (fT7)) . 
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GW frequencies MiXJ^m. = —'^(yim/^em^- For refer- 
ence also the orbital frequency MU is depicted on all 
panels. The vertical dashed line indicates the location of 
the maximum oi ft at t = tf2„ax ~ 4308. 4 A/, that corre- 
sponds to the time t = <lr, when the particle crosses the 
light ring tlr — 3M. Note again that we do not allow 
for any arbitrary time or phase shift between EOB and 
RWZ quantities. 

As previously observed poj the modulus of the in- 
splunge m = 2 waveform (left panel of Fig. is in very 
good agreement during the complete inspiral and during 
most of the plunge phase while overshooting the numer- 
ical one by about 10% around tLR. The frequency is 
indistinguishable by eye up to t « 4240 (relative error of 
0.5%), whereas it clearly underestimates the numerical 
one later on (by 11% at ^lr)- Similar results hold for 
the frequency of the m = 1 multipolc. In this case, how- 
ever, the insplunge amplitude is slightly smaller than the 
numerical one at <lr - 

As discussed in Sec. IIIIl a way to improve the in- 
splunge waveform in the strong-field-fast-motion regime 
is to include the NQC corrections, i.e. to consider the in- 
splumerg waveform, Eq. ([B]). The insplumerg waveform 
depends on the of™ parameters that must be tuned re- 
quiring some compatibility conditions with the numerical 
waveform. Previous work was mostly re- 

stricted to the quadrupolar m = £ = 2 waveform and 
considered only NQC amplitude corrections, i.e. fixed 
by construction in Eq. ([5]) (see, however. 



^22 



- ^22 



Ref. [10| for a preliminary investigation of the effect of 
for spinning binaries). The amplitude-related param- 
eters in /i^^ were generically fixed by imposing that (i) 
the maximum of the EOB waveform amplitude occurs at 
the time when the EOB orbital frequency peaks, 
and (ii) the maxima of the EOB and the numerical wave- 
forms amplitude agree at to„,ax- When only two param- 
eters are present [6[ , these conditions arc sufficient to fix 
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FIG. 3: Addition of NQC corrections and matching to QNMs; 1 = 2 multipoles. The (hght) dashed hnes refer to the bare 
insplunge waveform, without the addition of NQC corrections (dash-dotted line, blue online) nor of QNM ringdown (dashed 
line, red online). The vertical (light) dashed line indicates the location of the maximum of MQ.. 



with i = 1,2. Additional parameters are tuned by our conditions read 



locally fitting the numerical wavefori n M , fiot . Our pro- 
cedure builds upon previous works 0, [29l| with important 
differences. First, we take into account both phase cor- 
rections in Eq. The need for NQC phase corrections 
is motivated the frequency plots in Fig. [3] (note the rather 
large "gap" between the RWZ and the EOB insplunge 
frequencies at tLR) and requires two more conditions to 
fix the af™ parameters with z = 3, 4. Second, we deter- 
mine the four parameters af™ by demanding agreement 
of EOB and RWZ waveforms in amplitude, frequency, 
and their first derivatives at a given time t = t„i. Third, 
the procedure is applied to all multipoles. In formulas. 



EOB^+ ^ _ /I RWZ/, \ 



^EOB 



A^^^d \ — A 

{tm) = ■^m'^i^ra) 
(tm) — W^^^(tm) 
(tm) — l^frn'^i'^rn) 



EOB/ 



• EOB 



(14) 
(15) 
(16) 
(17) 



The necessary information extracted from the RWZ 
waveform is displayed in Table |IT1 The only remaining 
freedom in this procedure is the choice of t„i- Accord- 
ing to the usual EOB prescription we choose it from the 
EOB dynamics as the time ia„ax when the orbital fre- 
quency peaks, i.e. the time when the small black hole 
crosses the light ring at 3M. In our setting we have 
tm ^LR = io„,, = 4308. 4M. We note that does 
not exactly coincide with the time locations of the 
peaks of the Aim^s. In particular, A22 peaks 2.56M ear- 
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Her than tLR, while A21 peaks 9.40M later. Thus the 
maxima of the amphtude of the muhipolar waveform 
have no strict relation with the light-ring crossing, but 
they remain clearly identifiable points in the waveforms. 
An analysis of the amplitude maxima is reported in Ap- 
pendix |B] The information in Fig. [3] is completed by 
Table |V] in Appendix [BJ The third column of the table 
lists the time shift ^max^^LR with respect to the light-ring 
crossing iLR at which each multipolar amplitude peaks. 

The £ ~ 2 insplumerg amplitudes and frequencies 
with NQC parameters determined from Eqs. (|14p - PT)) 
are shown as dash-dotted lines (blue online) in Fig. |31 
The insplumerg is very effective in reproducing the nu- 
merical data. The corresponding maximum phase dif- 
ference accumulated between iLSO and amounts to 
0.025 rad for both multipoles, and the maximum rela- 
tive difference in amplitude is below 2% (see also Figs. |4] 
and [5] below) . 

The complete FOB amplitude and frequencies with the 
QNM contribution matched at <lr arc shown as dashed 
lines (red online) in the four panels of Fig. [3] (labeled as 
FOB: FuU). We match both muUipoles to three QNMs, 
using the values computed in [E^l- For m = 2 we use 
three positive frequency modes, 17^21 with n = (0,1,2) 
and we neglect the negative-frequency mode contribu- 
tion because it is very small [l^. For to = 1 we use two 
positive frequency modes, ct^jj with n = (0, 1) and one 
negative mode ctqj , to qualitatively reproduce the related 
oscillation in Muj2i. The center of the 3-point matching 
interval [t™ — A/2, im+A/2] ("comb") is chosen = ^lr 
and its width is A = M. The matching procedure is ro- 
bust for different choices of A and employs a minimum 
number of QNMs. A larger matching interval [2^ does 
not yield any improvement in the current setting, while 
a pointwise matching (A — >■ 0) leads to inaccurate re- 
sults. From the figure it is clear that the representation of 
the ringdown does capture the behavior of the numerical 
waveform more accurately in the to = 2 case and less ac- 
curately in the to = 1 case. This is a consequence of two 
facts. First, in the FOB framework the transition from 
merger to QNM ringdown is localized, by construction, 
at one single point t — tm- Second, we choose to deter- 
mine the NQC parameters and to match to QNMs at the 
same time, namely tm = ^lr- We have explored in Ap- 
pendix |2] the possibility of matching QNMs at a shifted 
time (common to all multipoles) tmatch > ^lr while keep- 
ing the tuning of NQC corrections fixed at ^merger = ^lr- 
The agreement of modulus, frequency, and phase during 
ringdown significantly improves for all multipoles, at the 
price, however, of one arbitrary shift parameter, namely 

At tmatch tmci-gcr- 



TABLE III: Values of the NQC coefficients af™ entering the 
h^2'^ factor, Eq. These numbers are obtained imposing 

the compatibility conditions Eqs. (|14|l - (|17| l between EOB and 
RWZ waveforms aX t = tm- 
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FIG. 4: Time evolution of the phase difference A0f ^^^^^^ = 
4)1^^ - '/'Sf ^ between the full EOB and RWZ multipolar 
waveforms. The dash-dotted vertical line locates the light 
ring. 



In the end we put together all the multipolar information 
to obtain the complete TZ{h^ — ihx) waveform given by 
Eq. ©. 



V. COMPLETE MULTIPOLAR WAVEFORM 

In this section we evaluate the performance of the pro- 
cedure discussed above to determine the NQC parameter 
fli when it is applied to the other multipoles with i > 2. 



A. Multipoles with £>2 

The qualitative behavior of the insplunge waveform 
during the late-plunge phase for ^ > 2 is analogous to 
the i = 2 case. The m = £ EOB waveform amplitude 
is always slightly larger than the corresponding numcr- 
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TABLE IV: Phase diflFerence 



BOB _ j,EOB J,RWZ 
= " Vim 

and relative amplitude diflFerence [AA/Af'^'^ = {Af°^ - 
Afm^)/Af^^ between the full EOB and RWZ waveforms 
computed at to ~ 500M and at the LSO crossing, tLSO = 
4076. lAf, for all the (m / 0) muhipoles with 2 < £ < 4. 
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FIG. 5; Time evolution of the relative amplitude difference 
(AAM)f™ = - Af^^)/ylf^z between the fuU EOB 

and RWZ multipolar waveforms. The dash-dotted vertical 
line locates the light ring. 



ical one around merger and progressively smaller when 
\ < m < Analogously, the gap between RWZ and 
EOB frequency is always smaller when m = i and pro- 
gressively larger when to decreases. This suggests that 
one should obtain numerical values for the a\^^ coeffi- 
cients that arc systematically larger when to decreases, 
as we actually find (sec Tabic IIIip . Note in addition that 
the value of is of order unity, while that of is of 
order 10^^ (as well as most of the others). These nu- 
merical values indicate the consistency of the procedure 
when it is applied to the EMRL case or to the compara- 
ble mass case. The following values were obtained in Q : 



J12 



-0.036 347 and af 

•21 



1.2468 in the equal- mass 
case {v = 1/4) and af ^ -0.017017 and af ^ 1.1906 in 
the 2:1 mass case {ly = 2/9). The variation of the coeffi- 
cients ap is mild when v varies between and 1/4, and 
it indicates that most of the v dependence is already ac- 
counted for by the linear momentum p^. and its first time 
derivative. It will be interesting to investigate whether 
this consistency remains (and at what accuracy) for the 
other multipoles. 

Figures 0] and [5] quantify the phase differences 
as well as the relative amplitude difference 
{AA/ A)i„i for all multipoles. These quantities remain 
quite small until the light ring: max(A0f °^^^^) < 0.15 
rad and |(AA/A)fo„| < 2.5%. Table |TV] (that is the ana- 
logue of Table |l| complements the late-time information 
given by Figs. S] and [5] by listing the numerical values of 
the phase difference and the relative amplitude difference 
at to = 500Af and at the LSO crossing. The accuracy ob- 



tained until the light ring is then lost during the ringdown 
part, especially for subdominant multipoles. Because we 
are determining the NQC phase-correction parameters 
(flg^jfll™) by means of one condition on tu = (j> and one 
on LOim = 4>irm thc phasc difference itself is not exactly 
zero at the matching point tm = ^lr (as is the case for 
the amplitude difference). This choice of imposing com- 
patibility on ugrn and (jj^m might look overcomplicated: 
One could just impose compatibility of the phase and 
its first time derivative. Working only with derivatives 
of the phase may allow us to extend our procedure to 
the comparable mass case using NR data to tune the 
NQC corrections. While in this work we have an unam- 
biguous correspondence between the dynamics and the 
numerical waveform, in the NR simulation this is not the 
case: The dynamics is not evidently available and one 
must typically rely only on waveform information that 
comes with some arbitrary initial phase. A procedure 
to fix (ag™, al™) independently of the GW phase is then 
preferable. In this respect Fig. [4] indicates that compat- 
ibility conditions on ugm and effectively yield phase 
differences that are quite small at the matching point 
(especially for the m = £ multipoles). As we shall see 
in the next section, the larger phase differences that are 
obtained for subdominant multipoles at [fikc the (4,1) 
case] have no practical influence on the complete wave- 
form. 

Finally we tested the procedure by using as the match- 
ing time thc location of thc maximum of the m = £ = 2 
amplitude t^^^, that occurs 2.56Af before tLR- We mea- 
sured thc corresponding values of the RWZ functions 
and computed thc NQC correction factors accordingly. 
Because t'^^^ ~ tLR, the EOB waveforms obtained in 
this way are fully compatible with those shown so far. 
Moreover, we also explored the effect of using a different 
matching time for each multipole, i.e. using as tm the 
location of thc maxima of the various Aim- For some 
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multipoles this procedure is very effective in obtaining 
an accurate representation of the modulus, hke, for ex- 
ample, in the (2,1) or (3,1) case, although it typically 
fails to reproduce the frequency. The reason for this is 
that in several situation the matching position already 
corresponds to the final growth of the frequency and the 
simple representation given by our NQC corrections does 
not capture the correct behavior there. Therefore, we 
prefer to use the simple procedure discussed so far, even 
if it might lead to a relatively inaccurate representation 
of the ring-down part of the amplitude for certain multi- 
poles (as for £ ^ 2, m = 1). 



B. The complete h+ — ihx waveform 

Now that we have assessed the quality of the EOB rep- 
resentation of the higher-order multipoles, let us compare 
the RWZ and complete EOB gravitational waveforms. 
This comparison is shown in Figs. [6] and [T] Figure [6] 
displays the late-time evolution of the two GW polariza- 
tions, TZh^/v (top) and TZh^/v (bottom), that are com- 
puted by summing together all the m ^ multipoles up 
to I = A using Eq. (jSj and considering a fiducial direction 
of emission (6*,$) = (7r/4,0). The corresponding phase 
difference is shown in Fig. [71 At the beginning of the in- 
spiral (top panel) the dephasing is as small as 1 x 10"'^ rad 
and is seen to progressively increase up to only 4 x 10^'^ 
rad at the LSO crossing (dashed vertical line). The phase 
difference increases then by an order of magnitude dur- 
ing the plunge phase and reaches about 0.03 rad at the 
light-ring crossing (dash-dotted vertical line). The main 
contribution to the phase difference comes from the i = 2 
modes and the m = £ modes with i > 2. The other modes 
have a much smaller impact (especially around merger), 
even in the EMRL. The relatively rough description of 
the ringdown structure (especially for the £ = 2, m = 1 
mode) implies that the phase difference is larger after the 
light-ring crossing, though it is mainly oscillating around 
zero. 



VI. CONCLUSIONS 

In this paper we have discussed the properties of 
the gravitational radiation emitted during the transition 
from quasicircular inspiral to plunge of two nonspinning 
black holes in the EMRL within the EOB framework. 
We considered for the first time the whole multipolar 
structure of EOB-resummed waveforms and we compared 
and calibrated them against recently calculated Regge- 
Wheeler-Zerilli numerical waveforms H^H^. The target 
numerical waveforms are extracted at null infinity via the 
hyperboloidal layer method [48l-l5l| . The binary dynam- 
ics is modeled in both cases for a point particle moving 
on a Schwarzschild background under the action of 0{v) 
dissipative radiation-reaction force computed using an- 
alytically resummed 5PN results [1, . As a paradig- 
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FIG. 6: Comparison between the RWZ and EOB com- 
plete gravitational waveforms taken in the direction (6, "l") = 
(0, 7r/4). The two panels show the time evolution of the two 
polarizations (/i+, hx) computed including up to ^ = 4 multi- 
poles for V = 10~^ (the m — Q multipoles are neglected). 



matic example, we have considered a binary with mass 
ratio V ~ 10"'^, initially at a relative separation of 7M, 
that inspirals for about 37 orbits before plunging into 
the black hole. The setup of our point-particle "labora- 
tory" is sufficiently general to allow us to gather infor- 
mation that can be useful both in the analytical models 
of waveforms emitted by EMRI (target sources for LISA) 
as well as for the coalescence of comparable mass black 
holes (target sources for ground-based interferometers). 
Our results can be summarized as follows. 

Quasiadiabatic inspiral. — At the beginning of the 
inspiral, the phase difference between the complete 
EOB and RWZ waveforms (computed without allow- 
ing for any relative time and phase shift) is very small: 
^^EOBRWZ _ 5 X 10-4 rad. This value is consistent 
with the estimated uncertainty related to the residual 
phases Sfrn entering the EOB waveform known only up 
to 4.5PN level. During the ^ 33 orbits of the in- 
spiral after the junk radiation (t < 500Af) up to the 
LSO crossing (corresponding to ~ 420 rad of total GW 
phase) the system accumulates only —2.5 x 10^'^ rad, 
i.e. A0EOBRWZ/0RWZ ^ _5 95 ^ iq-s. Such remarkable 

phase coherence that is obtained with the EOB insplunge 
waveform — without any tunable parameter — strongly 
indicates the aptitude of EOB waveforms to model EM- 
RIs for LISA. 

Our conclusions are compatible with those of Refs. [30l . 
[3l| although there are two important differences. First, 
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FIG. 7: Time evolution of the phase difference A^'^"''"^^^ = 
0BOB _^RWZ between the EOB and RWZ waveforms of Fig.[6l 
The top panel displays the behavior of this function from the 
very beginning (corresponding to initial separation r — 7M); 
the bottom panel focuses on the late-time behavior. 



our two waveforms are computed from the same dynam- 
ics in order to focus only on waveform comparison so to 
test the efficiency of the resummation of the EOB wave- 
form. Second, we do not further calibrate the i/ = 
EOB-resummed flux (and thus the radiation reaction T^) 
to circularized exact data [111- We believe that an ad- 
ditional tuning of higher PN contributions to the flux, 
though necessary for dynamical accuracy, would have 
only a marginal influence on our results. 

We finally remark that our setup and the accuracy of 
our data permitted us to assess the quality of the approx- 
imation to the (Sfm's residual phases entering the EOB 
waveform. Note that we have used the Sim^s in their 
standard Taylor-expanded form and we have not 

attempted to further resum them using nonpolynomial 
expressions. This might certainly be interesting to ex- 
plore to further reduce the (small) phase gap we have at 
the beginning of the evolution. 

Transition from inspiral to plunge, merger and 
ringdown. — For the first time we have explored the im- 
pact of NQC corrections to the complete multipolar 
waveform, including higher multipoles with 2 < £ < 4. 
The addition of NQC corrections is important to improve 
the EOB and RWZ modulus and phase agreement to- 
wards merger. We have proposed a simple procedure 
to determine NQC corrections on both the phase and 
amplitude for each multipole; four parameters are re- 
quired, two for the amplitude and two for the phase. 



They are determined by imposing compatibility between 
EOB and RWZ waveform amplitude, frequency, and their 
first derivatives at the light ring, i.e. the maximum of the 
orbital frequency. The procedure is robust and applies di- 
rectly to all multipoles (including those with £ > 4, that 
we have not explicitly discussed in the text). The com- 
plete EOB gravitational waveform (summed up to ^ = 4) 
shows a remarkably good phasing and amplitude agree- 
ment with the numerical one up to merger (±0.015 rad). 
After the light ring crossing we have a total phase dif- 
ference of about 0.25 due to the approximate treatment 
of the ringdown (via matching to QNMs), although it 
mainly oscillates around zero. The maximum relative 
amplitude difference is about 2.5% just before the light- 
ring. We emphasize that the exquisite phase agreement 
that we find at merger crucially relies on the calibration 
of NQC corrections to the phase. 

The procedure discussed here to determine the NQC 
corrections can be directly applied to EOB and NR com- 
parisons for comparable mass ratios generalizing current 
techniques. However, when p j^O, the procedure is more 
complicated due to the dependence of the EOB dynamics, 
notably of the Hamiltonian, on other flexibility param- 
eters that are also required to be determined (or con- 
strained) by NR data. In particular, one of the most 
evident physical effects entailed by these corrections on 
the dynamics is to displace the location of the "EOB- 
light ring" (i.e. the maximum of il) and thus the loca- 
tion of the matching time In our i/ = setting we 
can identify on the waveform unambiguously the time 
tLR that corresponds to the crossing of the light ring, 
because of the very good 10^'^ rad) phase ahgnment 
of the waveforms at early times and because the under- 
lying dynamics is the same. This allows us to measure 
the useful RWZ information at the correct location. As 
we emphasized in the text tm does not coincide with the 
time when A22 peaks, but it occurs 2.56M earlier. On 
the other hand, when dealing with NR data, the exact 
dynamics is not evidently available and one can rely only 
on waveform information. The peak of the exact orbital 
frequency, if it existed, should occur slightly after the 
peak of the A22 metric waveform amplitude. As a con- 
sequence, to apply the same discussed here to fix the 
NQC parameters and to keep the maximum of as the 
matching point, one should measure the four numbers 
per multipoles slightly (say by ^ IM) after the peak of 
A22- This method is different from current methods in 
EOB and NR comparisons, i.e. fix the NQC amplitude 
corrections imposing that the EOB and NR A22 peaks 
coincide at the maximum of the EOB orbital frequency. 
Even if this procedure is not a priori incorrect, we stress 
that if we were following this prescription in our setup 
we would have obtained a significantly larger phase dif- 
ference ( ~ — 0.2 rad), accumulated (starting from the 
10~^ level) in the last 50M before merger. This suggests 
that a more detailed analysis of the impact of NQC cor- 
rections to EOB waveforms in the comparable mass case 
might be needed in the future. 
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FIG. 8: Improvement of the EOB 1 = 2 ringdown waveform when matching the QNMs at fmatch = tLR + 3M. As in Fig. O 
the (light) dashed lines refer to the bare insplunge waveform, without the addition of NQC corrections (dash-dotted line, blue 
online) nor of QNM ringdown (dashed line, red online). The vertical dashed line locates the maximum of MQ.. 



As a last remark, we emphasize that extraction of nu- 
merical waveforms at null infinity convinced us to avoid 
any further (arbitrary) phase and time shift, providing 
clean information about the accuracy of the analytical 
modeling of the GW phase in the EOB waveform. On 
the contrary, we have shown that waveforms extracted 
at the finite radius = lOOOil/ yield initial phase differ- 
ences (with the EOB ones) during the inspiral that are 
~ 3 X 10~^ rad for the (2,2) multipole and about twice 
as much, ~ 5.5 x 10~^ rad, for the (2,1) multipole (and 
even larger for the subdominant multipoles [48|). This 
fact strongly indicates, once more, that in any EOB and 
NR comparison it is necessary to work cither with NR 
waveforms extrapolated to infinite extraction radius or 
evolved up to null infinity (sol [60l | . 
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Appendix A: QNM matching at fmatch > tLR 

In this appendix we briefly explore how the quality of 
the postmerger waveform depends on shifting the match- 
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FIG. 10: Maxima of the modulus of the RWZ master function 
for different multipoles. The dashed lines are obtained by 
fitting the data with Eq. (jBlJ. 



frequency modes. On the contrary, the amplitude differ- 
ences in the ringdown are still present. 

The influence that the match to QNMs at tmatch > ^lr 
has on the total waveform (summed up to € = 4) can 
be appreciated in Fig. [5] for the two GW polarizations 
(/i+,/ix) (top panels; contrast them with Fig. [5]) and for 
the total phase difference (bottom panel; contrast it with 
Fig. [71). Notably, in the latter plot the amplitude of the 
oscillation in the phase difference is more than a fac- 
tor of 2 smaller than the standard case. The result is 
mainly due to the relatively inaccurate representation of 
the ringdown of the (2,1) mode. 



ing to QNMs at tmatch > to,„ax ~ ^lR j maintaining the 
determination of the NQC parameters at t = ^lr- While 
the results of Sec. II VI and fVl are already satisfactory, we 
experimentally observe that displacing the matching to 
QNMs by At = tmatch — tLR ^ 3M produces a further 
improvement in the final waveform. This approach is 
also suggested by a careful inspection of Fig. |3] in the 
region around tLR,: It is clear that (for both multipoles) 
the EOB insplumcrg frequency and amplitude want to 
remain close to the numerical one even after t ~ tLR, so 
that we do not seem to be strictly obliged to just match 
the QNMs at t,„ = tLR,. The result of matching QNMs 
at a shifted time is shown in Fig. |8] for 1 = 2 (compare it 
with Fig. [3]). The improvement in the closeness between 
frequencies and amplitude during the ringdown is strik- 
ing especially for the m = 2 mode. An analogous result 
is obtained for the other multipoles, the more important 
improvement affecting the m = i modes. Concerning the 
£ = 2, m = 1 mode (and similarly for other odd-parity 
multipoles with i > 2) one succeeds in having a more ac- 
curate representation of the oscillation in the frequency 
due to the interference between positive- and negative- 



Appendix B: Multipolar hierarchy at merger 

In this appendix we study how the waveform ampli- 
tude peaks = ma,x{Aim) depend on the multipolar 
order (£, to). The numerical values of A'^^ present some 
regularities that it is worth discussing in some detail. In 
Fig. [TU] we show A^^^/i^ for all multipoles up to £ = 8. 
The corresponding numerical values are listed in Table W\ 
for completeness. The multipolar structure is such that 
the peak of the waveform amplitude is about 10 times 
smaller when the azimuthal number £ decreases by one 
(for a fixed to), while, for a given £, it increases with to. 
As an example, the maximum amplitude of the multipole 
(4,3) is larger than the (3,0) one. 

A fit of the Af^"" data of Table |V] with the functional 
form 

ln{AT^n = c^i£)m + C2i£)e (Bl) 

shows that the behavior of the A'^^ is approximately 
exponential in m for each £; for each value of m A'^^, 
scales approximately as e^^ when £ > 2; i.e. we have 
^im^ ~ e™~^^. The precise numerical values of ci{£) and 
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TABLE V: Information about the maxima of the modulus of the RWZ muhipoles. For each multipole we report the time shift 
between the time location of the maximum and the light-ring crossing time tLR, the maximum of the modulus of the RWZ 
waveform, its derivative at tLR, the frequency, and its derivative at tf^x- 



I 


m 


(iSSx — ^lr) 


A max 


^ max / j-^m \ 
^im V'maxJ 




Ma;fm(iSSx) 


2 





1.106x10^ 


0.527x10"^ 


-0.220x10"* 


0.790x10"^^ 


0.149x10"^* 


2 


1 


9.400x10" 


0.107x10"" 


0.258x10"^ 


0.291x10"" 


0.591x10"=' 


2 


2 


-2.560x10" 


0.295x10"" 


-0.260x10"^ 


0.272x10"" 


0.587x10"^ 


3 





1.036x10^ 


0.258x10"^ 


0.208x10"^ 


0.313x10"^" 


-0.739x10"^' 


3 


1 


1.054x10^ 


0.569x10"^ 


0.767x10"'' 


0.411x10"" 


0.679x10"^ 


3 


2 


6.840x10" 


0.182x10"^ 


-0.108x10"^ 


0.452x10"" 


0.159x10"^ 


3 


3 


1.000x10" 


0.515x10"^ 


-0.519x10"'^ 


0.453x10"" 


0.109x10"^ 


4 





1.508x10^ 


0.248x10"^ 


0.321x10"^ 


-0.523x10"^* 


0.190x10"^'' 


4 


1 


1.072x10^ 


0.488x10"^ 


0.572x10"^^ 


0.552x10"" 


0.371x10"^ 


4 


2 


9.520x10" 


0.165x10"^ 


0.158x10"'^ 


0.626x10"" 


0.196x10"^ 


4 


3 


7.200x10" 


0.497x10"^ 


0.686x10"'' 


0.637x10"" 


0.189x10"^ 


4 


4 


2.820x10" 


0.145x10"^ 


-0.225x10"'' 


0.634x10"" 


0.152x10"^ 


5 





1.424x10^ 


0.293x10"* 


-0.178x10"^^ 


-0.208x10"^^^ 


-0.244x10"" 


5 


1 


1.422x10^ 


0.562x10-* 


0.547x10"'' 


0.745x10"" 


-0.110x10"^ 


5 


2 


1.122x10^ 


0.186x10"^ 


0.390x10"'' 


0.807x10"" 


0.312x10"^ 


5 


3 


9.140x10" 


0.622x10"^ 


-0.501x10"^ 


0.801x10"" 


0.263x10"^ 


5 


4 


7.620x10" 


0.173x10"^ 


-0.302x10"^ 


0.822x10"" 


0.222x10"^ 


5 


5 


4.120x10" 


0.521x10"^ 


-0.268x10"*^ 


0.817x10"" 


0.189x10"^ 


6 





1.664x10^ 


0.388x10"^ 


-0.283x10"^ 


0.188x10"^'' 


0.299x10"^^ 


6 


1 


1.376x10^ 


0.758x10"^ 


-0.595x10"* 


0.880x10"" 


0.332x10"^ 


6 


2 


1.332x10^ 


0.246x10"* 


0.461x10"* 


1.021x10" 


0.102x10"^ 


6 


3 


1.128x10^ 


0.818x10"* 


-0.199x10"* 


1. 001x10" 


0.297x10"^ 


6 


4 


9.240x10" 


0.266x10"^ 


0.254x10"'' 


0.985x10"" 


0.288x10"^ 


6 


5 


8.080x10" 


0.697x10"^ 


-0.247x10"^ 


1.008x10" 


0.249x10"^ 


6 


6 


5.140x10" 


0.216x10"^ 


-0.550x10"^ 


1.002x10" 


0.222x10"^ 


7 





1.588x10^ 


0.598x10"'^ 


-0.332x10"* 


-0.929x10"^^ 


-0.256x10"^^ 


7 


1 


1.588x10^ 


0.112x10"^ 


-0.212x10"* 


1.063x10" 


-0.878x10"^ 


7 


2 


1.364x10^ 


0.365x10"^ 


-0.350x10"^^ 


1.186x10" 


0.416x10"^ 


7 


3 


1.252x10^ 


0.123x10"* 


-0.128x10"* 


1.189x10" 


0.330x10"^ 


7 


4 


1.128x10^ 


0.390x10"* 


-0.131x10"* 


1.186x10" 


0.319x10"^ 


7 


5 


9.460x10" 


0.124x10"^ 


0.505x10"* 


1.171x10" 


0.311x10"^ 


7 


6 


8.540x10" 


0.311x10"^ 


-0.238x10"^ 


1.195x10" 


0.273x10"^ 


7 


7 


6.000x10" 


0.993x10"^ 


-0.692x10"^ 


1.188x10" 


0.250x10"^ 


8 





1.746x10^ 


0.951x10"^ 


0.700x10"^ 


-0.774x10"^* 


0.223x10"" 


8 


1 


1.534x10^ 


0.181x10"*^ 


-0.336x10"^^ 


1.199x10" 


0.460x10"^ 


8 


2 


1.524x10^ 


0.579x10"'^ 


0.141x10"^ 


1.400x10" 


0.728x10"^ 


8 


3 


1.366x10^ 


0.195x10"^ 


0.205x10"^ 


1.386x10" 


0.364x10"^ 


8 


4 


1.242x10^ 


0.646x10"^ 


-0.699x10"^^ 


1.372x10" 


0.354x10"^ 


8 


5 


1.136x10^ 


0.197x10"* 


0.547x10"^" 


1.373x10" 


0.338x10"^ 


8 


6 


9.720x10" 


0.621x10"* 


0.396x10"* 


1.359x10" 


0.332x10"^ 


8 


7 


8.960x10" 


0.150x10"^ 


0.153x10"' 


1.383x10" 


0.294x10"^ 


8 


8 


6.720x10" 


0.490x10"^ 


0.238x10"' 


1.375x10" 


0.276x10"^ 
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TABLE VI: Fit of the maxima of the modulus of the RWZ 
for different i multipoles as a function of m. Tlie fit assumes 
the functional form Eq. (|B1|) . The 90% confidence interval of 
the estimate coefficients is reported as well as the coefficient 
of determination . 



i 




90% Conf. 


C2(£) 


90% Conf. 




2 


0.86 


[0.87 0.85] 


-1.497 


[-1.491 


-1.502] 


0.989 


3 


1.02 


[1.02 1.00] 


-2.016 


[-2.013 


-2.019] 


0.995 


4 


1.05 


[1.06 1.04] 


-2.116 


[-2.114 


-2.119] 


0.993 


5 


1.07 


[1.08 1.06] 


-2.128 


[-2.126 


-2.130] 


0.995 


6 


1.09 


[1.10 1.08] 


-2.112 


[-2.111 


-2.113] 


0.997 


7 


1.09 


[1.10 1.08] 


-2.082 


[-2.081 


-2.083] 


0.997 


8 


1.10 


[1.11 1.09] 


-2.051 


[-2.050 


-2.052] 


0.998 



C2 it) are listed in Table I VI I together with the coefheient 
of determination that measures the quality of the fit ^ . 
We finally note that the value for A^^^ jv that we 



obtain here is not dramatically different from the val- 
ues obtained when the masses m\ and mi are compa- 
rable. For example, the value {A^^'^ Iv)^^^) = 0.295 is 
about a 10% smaller than the corresponding numerical 
value computed the equal-mass case, q = mi /mi = 1, 
{Af^'')/v)q=i - 0.321 and only about 6% smaller when 
g = 4, {A^2 l^)q=i 0.31. As already pointed out in [Ij, 
one can fit the values of obtained from a few (ac- 

curate) numerical-relativity simulations as a function of 
77 = (1— 4i/), where v = mimi / (mi+mi)'^ is the symmet- 
ric mass ratio and reduces to = l^/M when mi mo. 
Assuming a linear behavior in 77, as it was done in [6|, 
one obtains A':^^^ ~ 0.321i/[l - 0.0899(1 - 4i/)]. How- 
ever, it turns out that a better representation of the data 
is given by a quadratic behavior in (1 — 4i/) of the form 
A™'^'^ ~ 0.321j/[1 - 0.0162(1 - 4i/) + 0.0792(1 - 4j/)2]. In 
the future it will be interesting to check the numerical 
accuracy of this relation versus other NR data, as well as 
to generalize it to the other multipoles. 
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